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Abstract 



We calculate dispersion relations and propagators for the a, w, p and S mesons in 
relativistic nuclear matter with a proton-neutron asymmetry. In addition to the a-m and 
5-p mixings already present in symmetric matter, mixing occur between all components 
of the a-co-S-p system when the proton fraction differs from 1/2. 



(N| ' 1 Introduction 

O" 

This work extends previous results on the dispersion relation of mesons in relativistic, dense, 
hot nuclear matter, to systems where there exists a proton-neutron asymmetry. We will 
be working in a model of the quantum hadrodynamics type (QHD) pj, where the nucleon- 
nucleon interaction is mediated by the exchange of u, a, ir ... etc. mesons. In this work, 
the dispersion relations of the mesons are calculated by performing a linear response analysis 
around an asymmetric Hartree gound state. The propagators are then obtained by inverting 
the dispersion relations. The resulting expression corresponds to the random phase approxi- 
mation. 

A new feature of the asymmetric matter is that new meson mixing channels open. A 
standard case of meson mixing is that occurring for the a-uj pair, where one meson is converted 
into the other by desintegration and recombination of a virtual NN loop. Moreover, the 
isoscalar mesons can couple with the isovector ones whenever the distribution function of 
proton and neutrons differ. This has been calculated earlier for the ui-p pair (see e.g. 0, H) 
in relation with the interpretation of dilepton production in heavy ion collisions M, 0] . The 
5-p mixing was recently investigated by Q in symmetric matter as a further mechanism 
contributing to dilepton production, and also by [J/J and ||. We present here a derivation 
and systematical study of the dispersion relation and propagators of the coupled neutral 
meson system u>-a-p°-5 . The charged mesons will be studied elsewhere. 

Our calculations may apply to the physics of neutron stars or of heavy ion collisions. 

Nuclear physics in terrestrial conditions is for the most part well described by the N = 
Z hypothesis. Nevertheless, new data begin to be available with sizeable deviations from 
symmetry, by performing experiments with nuclei at the border of the stability lines. Indeed, 
with the planned construction of radioactive beam facilities in several countries, exciting 
new possibilities will open of investigating the equation of state and dynamical properties of 
asymmetric nuclear matter. Such experiments will hopefully be able to shed some light on 
the behavior of the asymmetry energy at high density M. 



Among nuclear physics applications, one can think for example of observing the excitation 



of collective modes in neutron rich nuclei [1C]. Another application concerns the description 
of the isospin dynamics in heavy ion collisions at intermediate energies. In particular, one 
of the accessible experimental observables, the balance energy, defined as the beam energy 
at which transverse collective flow disappears, is valued as providing a clean way to estimate 
the influence of medium efects on the nucleon-nucleon cross section [11]. Previous work 



[12 1 has shown that, besides short range effects described the Brueckner approximation, one 
cannot neglect the screening of the nuclear interaction arising from the exchange of dressed 
mesons. An interesting phenomenon would occur if a zero sound mode is excited, since it 
would appear as a pole in the meson propagator and a corresponding resonant behavior of 
the cross section Jl^]. This analysis can be extended in order to include the description of 
isospin effects ||, |||, ||. 

Large neutron to proton ratio is the rule in astrophysical situations, where the matter 
is in /3 equilibrium. In cold neutron stars, the proton fraction is about 10 % or lower. In 
supernovae explosions and hot proneutron stars, the proton fraction is higher, about 30 to 
40 %, but deviations from the symmetric nuclear matter cannot be neglected. 

In the context of the physics of supernova collapse and the early stages of proto neu- 
tron star cooling, the screening of the nucleon-nucleon interaction has consequences on the 



neutrino opacities 17, 18, 19, 2C]. In dense and hot matter, the neutrino-nucleon scattering 
rate is modified by the correlations which affect the nucleon current. In the random phase 
approximation, the neutrino-nucleon differential scattering cross section is given by 

da _ g| E u , 1 T ( fr(R)^r n m 

dE v dn ~ 32^3 Ev 1 - e-* [VVww Lflu) {1) 

where 11^^ is the retarded polarization obeying the Dyson equation 

n w = U WW + ^WS D SSa^sw ( 2 ) 

W and S represent vertices for the weak and strong couplings respectively and Dss a p is 
the dressed meson propagator. This propagator involves the meson dispersion relation in the 
denominator. The magnitude of the scattering rate therefore depends crucially on whether 
or not zero-sound modes turns up in the meson propagation within the integration range. As 



an example, the Skyrme calculation of IS], which present a spin zero sound mode, predict 



an enhancement of neutrino opacities, whereas the relativistic approach of ][L7L 18, E(J, where 



this mode is absent, predict a reduction of the opacities. Further study is necessary. 

The calculation of [17, 18] did not take into account the full picture of meson propagation 
in asymmetric matter. We have applied [2C] the propagators derived in the present work to 
the determination of the neutrino-nucleon scattering rate, and find that the reduction effected 
by RPA correlations is less efficient when all mixing channels are taken into consideration. 

2 The Model 

2.1 Lagrangian density 

We will use a version of the quantum hadrodynamics model developped by Walecka and 
coworkers in which the nucleon-nucleon interaction is mediated the six mesons taken into ac- 
count by the Bonn potential model, i.e a, to, tt, p, n, rj and 5 mesons. The latter two mesons 
are often neglected in mean-field type calculations by arguing that they lead to comparatively 



small corrections to the nucleon- nucleoli interaction potential in the vacuum [21, [2^]. Nev- 
ertheless, the 5 meson can play a significant role in the context of dense asymmetric matter, 
since the value of the 5 mean field controls the mass difference between proton and neutron. 
This field modifies the behaviour of the neutron star equation of state |)23|| . The value of 
the 5NN coupling has a strong influence on the proton fraction of matter in (3 equilibrium. 
Another justification is the fact that an extrapolation of Dirac-Brueckner-Hartree-Fock cal- 
culations to asymmetric matter was shown to be equivalent to a density-dependent mean 
field theory with a delta meson p4|, 25, |2q|. 



The 7T and r\ will mix with each other, and separate from the uo-a-p-5 sector. We will not 
study the tt-tj system here, since our main purpose was in fact to apply the results derived 
in this work to the calculation of the RPA reduction to the neutrino-nucleon scattering rate, 
where the pion does not contribute |2^] . The calculation of the mixed ~K-r\ dispersion relations 
would not anyway present any particular difficulty, using the same method as applied here 
to derive the equations which describe the a-5 subsystem. 

We start from the Lagrangian density 
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The hat on the Wigner operator and meson fields means that they are operators. The pion 
does not appear in this Lagrangian because it will not be investigated further in this work. 
Terms such as a-fr 2 or p .(d^if x 7r) also do not appear, beacause we decided to focus in this 
work on meson mixing originating from NN loops, and will not treat meson loops. A more 
realistic treatment should include them as well. 

This Lagrangian is non renormalizable stricto sensu due to the derivative couplings of 
the p meson (and ir if included with a pseudovector coupling). Nevertheless, this coupling 
is necessary for a good description of experimental data. In fact, renormalizability is not a 
necessary requirement for effective theories |27j]. Anyway, a renormalization procedure should 
be applied at the given level of approximation in order to avoid pathological behavior of the 
dispersion relation. To this end, we will add counterterms to the Lagrangian. A discussion 



of this issue is given in section §3.2 and Appendix B. 



2.2 Relativistic Hartree equilibrium 



The approach taken by [28| is applied here. The main steps leading to the dispersion relations 



are only briefly recorded; the reader is referred for details to [28|. It consists of writing a 
kinetic equation for the Wigner operator, and perform a linear response analysis of pertur- 
bations of the system around the equilibrium as obtained in the Hartree approximation. We 
first define the Wigner function 
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In the Hartree approximation, we suppose that the meson field operators a, to ... can be 
replaced by their classical expectation values an, ton ■■■ and that the equilibrium is homo- 
geneous. We will moreover assume that there does not exist charged condensates (only the 



third component of 5h and pn are non zero), so that the proton and neutron states remain 
well defined and the Wigner function is diagonal in isospin space. The nucleon fields are then 
plane wave solutions of a Dirac equation with effective momenta and masses. The Wigner 
function is obtained by introducing these solutions in the definition and taking the statistical 
average. It is given by 
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The middle two lines define the effective masses and momenta of the neutron and proton in 
terms of the classical fields. They are self consistent relations which can be solved for a given 
total density, proton/neutron asymmetry and temperature. 

The equation of state is obtained by calculating the energy momentum tensor and con- 
served baryon current 
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For further details on the equation of state, the reader may wish to consult the work of Kubis 
and Kutschera [p3| . 

The proton fraction is defined as 

n. 



Y P = —f— (13) 

Tin t Tip 

It can be given as an input parameter or be determined from the (3 equilibrium condition. 
The first case would be relevant for heavy ion collisions, where the asymmetry is determined 
by the choice of the nuclear isotopes participating to the reaction. In neutron stars and 
supernovae on the other hand, the timescales are long enough for the chemical equilibrium 
to establish through inverse (3 decay 



p + e — ► n + v 
yielding a condition on the chemical potentials which determines Y p 

M — P"n Mp — Me Mf 



(14) 
(15) 



In cold neutron stars, the neutrinos escape freely from the system, and we can moreover set the 
chemical potential of the neutrinos to zero. The proton fraction is of the order of 0. 1 depending 
on the density which is considered. In the last stages of supernova collapse and early phase 
of protoneutron star cooling on the other hand, the matter is opaque to neutrinos and \i v 
is determined by dynamical equations governing the lepton fraction Y^ = Y v + Y e . Charge 
conservation imposes Y e = Y p . For a typical value 1^=0.4, we obtain \x v = {^tv^hbYv) of 
the order of 200 - 250 MeV and Y p ~ 0.3 - 0.36. 

An example of the thermodynamics which are obtained from this model is shown on Fig. 
1 , where the effective masses and proton fraction are represented as a function of density for 
neutrino-free matter in beta equilibrium. It was calculated with parameter set fl49|) which 
fits the saturation point. Due to the coupling to the 5 field, the proton effective mass is 
higher than that of the neutron in this model. For the same value of the asymmetry energy 
at saturation, larger proton fractions are obtained at high density with a finite value of the 
coupling to the delta field, e.g. g$ = 5, than g$ = 0. 



3 Dispersion relations 

3.1 Linear response 

The Hartree solution was a particular approximation to the general evolution equations 

.7-9 
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The hat on the Wigner operator and meson fields means that they are operators. The 
physical quantities are extracted by taking statistical averages. Up to now, these equations 
are exact. A first approximation will be made by neglecting correlations, i.e. replacing the 
statistical average of products of meson operators with the Wigner operator, by the product 
of statistical averages: < <&F > — > < <I> >< F >. 

Let us now consider a small perturbation around the Hartree equilibrium. We will develop 
< F(x,p) >= Fh{p) + Fi(x,p), < a(x) >= an + o'i(x) ... etc. The Hartree equilibrium 
is obtained at zeroth order. At first order, we obtain equations for the perturbation which, 
after performing a Fourier transformation, are easily solved. The dispersion relations of the 
mesons are obtained by replacing the solution for the perturbation to the Wigner function 
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in the linearized, Fourier-transformed field equations for the mesons. For example, the dis- 
persion relation for the sigma meson reads 



(-q 2 + m 2 + 2bm<JH + 3caj I ) o\ 
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This defines the polarizations. The full set of dispersion relations is given in matricial form 
in § |3.3| . The index a is an isospin index. When taking the traces, it can be seen that the 
dispersion relations of the neutral mesons a, u>, p°, 5°, tt° decouple from those of the charged 
mesons p^ , 5^ , n^ . In this paper we limit ourselves to the dispersion relations of the neutral 
mesons; the case of the charged mesons will be the subject of future work. 

The nonlinear a 3 , a 4 couplings contribute terms on the left hand side which we can absorb 
in a definition 

M 2 =ml + 2bma H + 3ca 2 H . (26) 

We point out again that there are no meson loops (such as would be contributed e.g. by a 
7T-7T loop insertion in the a propagator) at this level of approximation in our formalism. They 
were discarded when we made the assumption < &F > < $ >< F >. One could obtain 
them by restoring meson-meson correlations, or insert such loops by hand basing ourselves 
on the diagrammatic approach . This will not be done in this work, since our first aim was 
to study meson mixing originating in NN loops. 



3.2 Polarizations 

Let us continue on the example of the a meson. Explicitely, the il^o- polarization is given by 

n = n^ + n( n ) 
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At zero temperature, the result obtained through the application of the linear response anal- 
ysis coincides with the one-loop approximation. The other polarizations are given by similar 
equations. Taking the traces and integrating over the angles, we finally arrive at the for- 
mulaes given in the Apendix. Some remarks are in order about the imaginary part at finite 
temperature and about the contribution of vacuum fluctuations. 



We will have to choose a prescription to go around the poles of the propagator Eq. (I 



In [28], the prescription was chosen as follows 
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with <pi(p) = j^yS(P 2 - M?)fi(p) 

At finite temperature, the real time Green's function formalism |29| defines several polar- 
izations depending on the position of the time arguments on the Keldysh contour. Their 
real parts coincide, the imaginary parts are related to each other by relations such as 
Im Tl R (w, k) = tanh (^r) 1m II 11 (w, k). It can be verified that the usual retarded polar- 
izations (2^] can be obtained by taking ±i e sign(po) instead of ± i e in the above expression 
(p7[). The imaginary parts given in the Appendix as well as those used in the calculation of 
the propagators will be those of the retarded polarizations. 



The last term 6(—po) in the expression of the Wigner function Eq. (10) describes the 
contribution of the vacuum. When taken under the integration sign in the expression of the 
polarizations, it yields a divergent term. Some renormalization procedure has to be applied 
in order to remove the divergence. It is not sufficient to ignore this term (by performing a 
normal ordering), because the subtracted term depends on the effective masses and therefore 
on the thermodynamical conditions. In the following we will apply dimensional regularization 



and counterterm subtraction as described e.g. in |3C, 25, |3T]]. What is new with respect to 
these references is that, due to proton-neutron asymmetry, the contributions of the mesons 
cannot be renormalized independently. 

The mixed polarizations IT^ and 11^ are both proportional to rf" = u tJ- — (q.u/q 2 )q^. They 
must vanish in vacuum, since the hydrodynamical velocity quadrivector is undefined in this 
case. Therefore, we can treat the renormalization of the a-6 subsystem independently from 
that of the uj-p subsystem. We review in this section the main steps of the renormalization 
of the a-5 subsystem and give the expressions of all renormalized vacuum contributions in 
Appendix B. 

After performing the dimensional regularization, the divergent contributions are extracted 
in the form (M 2 + M 2 )/e and q 2 /e in H"™ c or W 5 f, and in the form (M 2 - M 2 )/e in IT^ 
or IT^ C . The fact that M 2 + M 2 = 2(m - g a a) 2 + 2g 2 5 2 and M 2 - M 2 = 2(m - g a cr)gs5 
dictates the form of the counterterm Lagrangian one has to introduce, in order to neutralize 
these divergences. 

£ct = Z a d^o + l -A x a 2 + X -A 2 a 3 + ^A 3 a 4 

+Z 5 d^Sd^d + B 1 5 + B 2 ab 2 + B 3 a 2 5 2 + C 5 4 (28) 

Each of the couplings of the counterterm Lagrangian (Zi, Ai, Bi, C) is split (A\ = infinite 
+a%, etc) in an infinite part, which is chosen so as to cancel the divergence, and a residual 
finite part (zi, aj, bi, c), which has to be determined by imposing some physical conditions 
that the renormalized polarizations should verify. 

There are several possible choices, each of which defines a renormalization procedure. 
One possibility is to impose that the polarization and its first derivatives with respect to a 
and 6 vanish at some point, for example on the mass shell q 2 = m 2 (as e.g. in the work of 
Diaz Alonso et al. |28|]) or at q 2 = (as in the work of Kurasawa and Suzuki [p2| ). Another 
possibility is to require that the initial structure of the expression in terms of the effective 



masses is preserved (The procedures of e.g. Hatsuda et al. |33[ or Sarkar et al || do have 



this property, as well as "scheme 3" of |H}] .) Here we will work with the scheme of Kurasawa 
and Suzuki (called "scheme A" in the following) and scheme 3 of 1 31 ] called here "scheme 
B"). 



As was already pointed out in [31|, the final result, and in particular the behavior of the 



effective meson masses as a function of density, is very sensitive to the choice of renormaliza- 



tion procedure. We must consider this as an unsolved problem for the moment. There is some 



hope that one could solve it by applying "naturalness" arguments [27] and considerations on 
the symmetries of the underlying more fundamental theory, of which the present effective one 
represents a low-energy approximation. 

Actually, the asymmetric case offers us an example of how symmetry arguments do im- 
pose some further restrictions on the choice of the subtraction procedure. As a matter of 
fact, the mixing between the a and S mesons requires that the same counterterms should 
simulaneously cancel divergences in various polarizations. For example, the B3 is used to 
regulate the 5 2 /e divergence in H a a, the a 2 /e divergence in II^ and the aS/e divergences 
in n ff 5 and Ilsa- The B2 regulates divergences S/e in Ii a s and Hs a , and a/e in II^. The 
compatibility will impose additional relations between the finite 627 b^, etc. constants. It 
was found that the renormalization scheme used by |2^] must be discarded on this basis, 
since the additional constraints necessary for the compatibility of the counterterms are not 
fulfilled. Both scheme A and B used in this paper fulfill the compatibility relations. Scheme 
B preserves the structure of the expression as a function of the effective masses M p , M n but 
gives after renormalization n CT( $ 7^ n^, whereas these two quantities were given by the same 
expression before renormalization. For scheme A the reverse occurs, we have II^ = li^a after 
renormalization but the structure in terms of M p , M n is not maintained. 

Further details are given in Appendix B. 

3.3 Dispersion relations 

The equations for the dispersion relations can be cast in matricial form: 
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In a referential where the momentum transfer is equal to q^ 
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(uj, 0, 0, k), the T> matrix has 
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with the following definitions of a, b, c, d ... etc. 
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In writing the preceding equations, we used the decomposition of the polarizations on a basis 
of orthogonal tensors 
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If the non linear couplings b, c are non zero, one should replace everywhere m 2 by M 2 as 
defined in Eq. (J26|). Taking the determinant of this matrix, we obtain the dispersion relation. 



It is of the form det(T>) = m m 2 8j, 5l = 0, the transverse propagation modes being given 
by 



5 T = (q 2 - m 2 + U pT )(q 2 - m 2 + n wT ) - II 2 ^ 
and the longitudinal mode by 
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4 Propagators 



The propagators are obtained by inverting the dispersion relation Q = T> 1 . The propagator 
matrix has the same structure as the dispersion matrix. 
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The scalars isoscalar G a and isovector Gs propagators are given by 
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Mixed a-8 meson propagation in asymmetric matter is given by 

G a 8 = t- {-{q 2 -ml +Ii ujL )IispU ap rj 2 - (q 2 - m 2 + Il pL )Us UJ Il (7LU i] 2 



(32) 



(33) 



+ {U-su n CTp + lisp n CTW ) n pwL if 

- (q 2 ~ml+ Ii wL ){q 2 -m 2 + U pL ) - U 2 pujL \ U aS ] 



(34) 



We have the usual cr-cu mixing in the isoscalar sector and 5-p mixing in the isovector sectors. 
Morover, in asymmetric matter there are also mixing between the isoscalar and isovector 
mesons 5-uj and a-p 






if 1 



{(U 5u) U ap - U Sp U ffU! ) U Sp if - (q 2 -m 2 p + U pL )(q 2 - m\ + 11*) U aiV 

(q 2 -m 2 + n pL ) n fc n CT5 + (q 2 -m 2 + lis) n CTp u pu)L - u Sp n aS iv,l} (35) 



G S P = -y^s^p-JlspR^Il^Tf-tf-ml + U LjL )(q 2 -ml + U a )Us p 

+(q 2 -ml + n CT ) u Soj Hp^L + (q 2 - ml + n^) n CT< 5 n CTp - u aS n CTW n pwL }(36) 

G£ p = -^-{-{U SLU U (7p -ns p U (7W )Us iuV 2 -{q 2 -ml+U ujL )(q 2 -m 2 s+Us)U (Tp 

+(q 2 ~ m l + Hujl) n,5 P u aS + {q 2 -m 2 + u 5 ) n CTW n pwL - u Suj u a5 n pwL j (37) 

G&, = -^{-(U 5 ^U ap -U Sp U^)U ffp r I 2 -{q 2 -m 2 p + U pL )(q 2 -rnl+U tT )U Suj 

+ {q 2 -m 2 a + U a ) U Sp U pu)L + (q 2 - w 2 + U pL ) U aS U aiV - U ap U aS n pwL } (38) 
The propagators of the vector mesons uj and p can be decomposed as 



G^ = GulL^ + GutT^ + GuqQ^ 



G% 



G p l L^ v — G p t T^ u + G p q Q 
GoujL L^ + G 0(jJ t T^ 



I lis 



(39) 
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with for the u meson 

a 



— {(q 2 -m 2 5 +n 5 ) U 2 ap i] 2 + (q 2 - m 2 a + II CT ) U 2 Sp V 2 - 2 Ii Sp U ap Yl a5 tf 
+ (q 2 -rf, + n pi ) [{q 2 -m 2 a + U a )(q 2 - m 2 s + 11*) - Ii 2 5 ] } 



GujQ 



- — (q 2 -m 2 p + U pT ) 



1 



mf 



(40) 
(41) 

(42) 



for the rho meson 

G pL = ^{(<i 2 - m2 + u s) n L v 2 + (q 2 -ml + n a ) ni r? 2 - 2 n^ u auj n ff5 v 2 

Hq 2 ~ml + U wL ) [(q 2 -m 2 a + U a )(q 2 - mj + U 5 ) - n 2 aS ] } 

- — (q 2 -ml+Ii ulT ) 
Ot 

1 



G p T 
GpQ 



111- 



(43) 
(44) 

(45) 



and finally for the mixing between the to and p mesons 

Gw, - —\-{q 2 -ml + Il a )Tl Sp n 5uJ i] 2 - (q 2 - m 2 + Ii 5 )Ii ap 11^ r? 2 



T pLxjL 



Sl 



+(U 5uJ U ap + Hsp u aul ) u a5 v 2 

- \q 2 -ml + U a )(q 2 - mj + 11*) - 11^] U^l] 



G 



n 



pu)T 



puT 



St 



(46) 
(47) 



5 Numerical results 
5.1 Choice of parameters 



At the mean field level, the equation of state is determined by equations (11,12). The cou- 
pling constants were adjusted in order to reproduce the properties of the saturation point. It 
is enough to adjust five parameters in order to reproduce the five basic experimentally mea- 
sured properties: saturation density n sat , binding energy B/A, incompressibility modulus K, 
effective nucleon mass m* and asymmetry energy a a- Usually these parameters are taken to 
be the meson- nucleon coupling constants g a , g^, g p and the non linear sigma couplings b, c; 
and the coupling of the 5 meson is set equal to zero. The meson masses are kept equal to 
their value as quoted in the Particle Data Book (m w = 782.6 MeV, m p = 769 MeV, m.^983 
MeV), and the a meson mass is taken to be equal to the standard result of the Bonn potential 
model m a = 550 MeV. The tensor coupling of the rho meson is derivative and does not enter 
in the equation of state at the mean field level. It will be taken equal to the Bonn potential 
value f p /g p = 6.1 or to the vector dominance model value f p /g p = 3.7. 

If the S meson coupling is not equal to zero, we have one more parameter to play with. 
When keeping the parameters g a , guj,b, c constant, varying g$ rescales the asymmetry energy, 
which has to be readjusted by changing the value of g p . The 5 is usually discarded, moreover 
some studies concluded that is does not noticeably affect the description of stable nuclei [|34j . 
On the other hand, interest in the 5 meson has been revived since it appears that density 
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dependent field theory needs a 5 with a large value of the coupling in order to reproduce 
Dirac-Brueckner calculations of asymmetric matter [|3, B6|. Whether or not a <5 meson is 
present can be important for the equation of state of neutron stars (£3|. Among other things, 
a nonzero coupling to the 5 manifests itself by a splitting of the proton and neutron effective 
masses. Hopefully, experiments planned with neutron rich nuclei will help clarifying this 
issue. 

In principle, renormalization of vacuum fluctuations appearing in the expressions of the 
nucleon effective mass and energy density should be performed, yielding the relativistic 
Hartree approximation. The procedure is standard and does not present difficulties (see 
e.g. p8f| ). However it is known [|], [35|] that, after effecting a readjustment of the coupling 
constants, the equation of state is very similar in the renormalized (Hartree) and unrenor- 
malized (mean field) case. Since the thermodynamics enters the polarizations only through 
the effective masses M n , M p and chemical potentials fi n , fj, p , we take the pragmatic point of 
view that it will not change the dispersion relations whether the underlying thermodynamics 
were generated by the mean field or Hartree approximations, and limit ourselves to fit the 
mean field. Reasonable values for the experimental data are 

n sat = 0.17 fm" 3 , B/A = -16 MeV, K = 250 MeV, m*/m = 0.8, a A = 30 MeV (48) 

and can be reproduced with the following parameter set 

^ = 8.00, 9uJ = 7.667, b/gl = 9.637 10~ 3 , c/g 4 a = 7.847 HT 3 

(g p ,g s ) = (3.685,0.) or (5.203,5.) (49) 

Note that the parameter c is positive, so that we are safe from any unpleasant instabilities 
which may appear |j(J in commonly used parametrizations where this is not the case. 

This parameter set will be referred to in the following as set (1). For the dispersion 
relations we will use moreover other parameter sets besides this one. As a matter of fact, at 
this level of approximation, there is no compelling reason why the couplings describing the 
thermodynamics should be the same as those entering the RPA. When the dressed meson 
propagator is used in calculations of the RPA corrections to neutrino opacities, the scattering 
processes are usually described by parametrizations such as given in Bonn potential model. 
For example, from [22| 

= 10.20, g u = 15.85, b = 0, c = 

(50) 
A w = 1500 MeV 

The cutoffs A CT , A w , A p and A5 enter the definition of form factors which multiply the coupling 
constants in the calculation of the polarizations. They are introduced in order to take into 
account the effects of the finite size of the nucleon. 

/A 2 -m 2 Y a 
g a ^ g a ( ," 2 _ £ I , a£ {a,uj,p,5} , n a = l except for n p = 1 or 2. 

This parameter set will be referred to in the following as set (2). 

From ]37j we can take the parameter set of Brockmann-Machleidt potential C which also 
fits the free nucleon-nucleon scattering data, but with a higher value of the coupling to the 5 
meson 

g a = 10.044, g w = 15.853, 6 = 0, c = 

g p = 3.455, f p /g p = 6.1, gs = 7.985 (51) 

A CT = 1800 MeV, A w = A s = 1500 MeV, A p = 1300 MeV, n p = 1 
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9a " 


= 10.20 


9lu = 


15.85, b 


= 0, c 


9p~- 


= 3.19, 


fp/9 P 


= 6.1, gs 


= 3.73 


A a 


= A 5 = 


A p = 2000 MeV, 


n p = 2 



This parameter set will be referred to in the following as set (3). 

A way out of the dilemna between chosing parameters which fit the saturation point, or 
the scattering data in vacuum, could be working with one of the recent parametrizations of 
the results of Dirac-Brueckner calculations (which obtain the saturation point using the free 
potential as input) in terms of density-dependent coupling constants. A promising model 
is the density dependent relativistic hadron field theory (DDRH) which is lately beeing de- 



veloped by Lenske and collaborators [25|. A serious study would involve performing linear 
response analysis on the field equations obtained from the DDRH lagrangian density. This 
is currently under consideration. 

5.2 Transverse modes 

Depending on the values of the coupling constants and cutoffs on the one hand, and of the 
renormalization scheme on the other hand, we may have a different structure of the branches. 
If the renormalization is not performed, we have spurious branches, as is known from previous 
studies [^, [H|]. In particular, when we define the effective meson masses as the solution u 
of the dispersion relation on the axis k = 0, a kink structure is found, whatever strong form 
factor is applied. When the renormalization is performed, the kink disappears and normal 
branches emerge. 

In asymmetric matter, the transverse mode of the omega meson mixes with that of the 
rho meson. With renormalization scheme A (using the method of Kurasawa-Suzuki), we 
obtain no more branches than the normal ones with parameters sets (2) or (3). On the 
other hand, with parameter set (1), there is also a zero sound branch at finite momentum, 
similar to the one conducing to the (spurious) pion condensation. This zero-sound mode is 
related to the high value of the rho meson coupling constants necessary in order to reproduce 
the correct asymmetry energy, together with the high value of f p /g p = 6.1 of the Bonn 
potential. As a matter of fact, the zero-sound mode is appreciably reduced when we take 
the vector dominance model value f p /g p = 3.7. It disappears completely for f p = 0. On the 



other hand, recent calculations from QCD sum rules [38] would also favor a high value of 
fp (fp/dp = 8.0 ± 2.0). The zero-sound mode is also present with parameter set (1) when 
using the other renormalization scheme (B), although much weaker. When present, the zero 
sound mode is somewhat reduced at larger proton-neutron asymmetry. In any case, it will be 
possible to eliminate this spacelike mode with a contact interaction |8| of the Landau-Migdal 
type (as in the ir-p +g' model), which also helps to correct the singular 5(r) behavior of the 
rho exchange potential at short distance. 

The expression of the omega polarization is identical in schemes A and B. On the other 
hand, the result for the rho meson is quite different depending on the chosen renormalization 
scheme. As shown is pl| , the effective p mass increases with increasing density in the renor- 
malization scheme of Kurasawa-Suzuki (A), whereas it decreases in renormalization scheme 
(B) which keeps the formal strucutre of the vacuum terms as a function of the effective 
nucleon masses M n , M p . 

The position of the normal branches changes with asymmetry. It is best seen by plotting 
the variation of the effective rho and omega meson masses with asymmetry. The effective 
meson masses are defined as the solutions uj of the dispersion relation 8t{w, k) = at vanishing 
three momentum k. They are shown on Figs. 2 and 3, (rightmost upper and lower panels), 
at four times saturation density and vanishing temperature, with parameter set (1). The 
continuous lines were obtained by keeping the 10-p mixing while the dashed line is obtained 
by setting it to zero. 

In renormalization scheme A, the omega meson effective mass is the lowest of the two 
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and decreases with increasing asymmetry (= smaller Y p values), whereas the rho meson mass 
increases with increasing asymmetry. In renormalization scheme B, the reverse occurs. The 
effective rho mass is now the lowest and decreases with increasing asymmetry. As in the 
case of the behavior of the p mass with density pTJ, the uncertainty in the choice of a 
renormalization procedure once more precludes a reliable prediction. 

In order to evaluate the strength of the mixing, we can define the mixing angle as in 0] 



tan(2 9 wpT ) = 2 ^^ T (52) 

m z p -m^- Ii pT + II^t 

The dependence of Q^ p t in the density and asymmetry parameters is represented on Fig. 4. 
The curves are labelled by the value of the proton fraction Y p . The mixing angle is calculated 
at a fixed three-momentum k = 500 MeV, at the values of uj which are solutions of the 
dispersion relation. We have two angles, one for each meson branch. The mixing angle Q^pT 
vanishes in symmetric matter and/or at zero momentum. Non negligible values are obtained 
in asymmetric matter, especially if renormalization scheme A is used. 

5.3 Longitudinal modes 

When studying the solutions of the longitudinal dispersion relation 5l, we obtain the normal 
branches corresponding to each of the four mesons involved in the longitudinal part of the 
dispersion relation. Again, the behavior of the corresponding effective masses as a function of 
density depends on the choice of the coupling constants, cutoffs and renormalization scheme. 
In both schemes, the sigma meson mass increases with density, more rapidly in scheme B. 
The rho meson mass increases in scheme A and decreases in scheme B. The delta meson mass 
decreases in scheme A and increases in scheme B. The omega meson mass coincides in both 
schemes, it first decreases with density and then starts increasing. 

The behavior of the effective masses at a fixed density ub = 4 n sa t as a function of the 
proton fraction Y p is represented in Figs. 2 and 3. The effective masses obtained for the 
omega and rho meson from the longitudinal dispersion relation coincide with those obtained 
from the transverse modes, as it should be. The delta mass is larger at smaller Y p in both 
schemes. The sigma mass increases in scheme B but first decreases and then slightly increases 
with asymmetry in scheme A. 

Besides the normal branches, we have also heavy rho meson branches (see also e.g. pl|) 
which would disappear if a stronger cutoff were applied. These will not be further studied 
here. Finally, there are zero sound branches. The zero-sound branch due to a-uj mixing is 
already well known from studies of symmetric matter. It is present in both renormalization 
schemes with parameter sets (2) and (3). 

In Fig. 5 we show the dependence of the zero-sound mode due to a-uo mixing with the 
asymmetry, at four times saturation density and vanishing temperature. It is seen that at 
high asymmetry the strength of the zero sound is only slightly reduced. The associated 
characteristic velocity increases, resulting in a small rotation of the zero sound branch in the 
u>-k plane. When studying the imaginary parts of the polarization, it is seen that at zero 
temperature the lower branch of the zero sound mode is suppressed by Landau damping, 
whereas the higher branch is not. 

In the delta-rho sector, no zero sound branch similar to that existing in the o-io system 
was found, whatever parameter set was used. Some spurious branches may appear in the 
spacelike region at high momentum, if we work with parameter set (1) and in renormalization 
scheme A. They disappear when applying a moderate cutoff. 
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5.4 Effect of the mixing on the propagators 

In the preceding paragraphs, we obtained results for the effective meson masses and mix- 
ing angles, which concern the timelike region of the u-k plane. In order to characterize the 
influence of meson mixing in the spacelike part of the u-k plane, we also studied the prop- 
agators at fixed k for varing u with the condition to < k. The knowledge of the behavior 
of the propagators for this parameter range can be of use for the study of RPA corrections 
to neutrino-nucleon scattering, since it enters the definition of the RPA correction to the 
polarization insertion of Eq. (||) . 

A sample of the results is displayed in Figs. 6 to 8. The thermodynamical conditions 
were chosen to be tib = 4n sat , T=20 MeV, neutrino free matter in (3 equilibrium Y v = 0. 
The exchanged 3-momentum was fixed at k = 50 MeV. Similar features were observed for 
other values of ns and k. 

The real and imaginary part of the a meson propagator are shown on Fig. 6. The G auJ 
and G u l behave in the same way. The behavior of the real part of G a is dominated by the 
pole at the location of the zero sound branch. The lower part of the zero sound branch falls 
into the region of Landau damping (u/k = 0.67 on the figure) so that it does not give rise 
to a pole. On the other hand, the upper branch falls on the fringe of the Landau damping 
zone at zero temperature, so that the corresponding pole {u/k = 0.83 on the figure) remains 
largely unscreened. The effect of mixing is not very strong. 

As we saw in the preceding paragraph, there does not exist a 5-p zero sound branch. Due 
to the full a-u-p-5 mixing however, the propagators of the 5-p subsystem get contaminated 
by the zero sound pole in the a-u subsystem, as can be seen on Fig. 7, where we plotted the 
Gs p propagator. The G$ and G p l propagators behave in a similar way. 

Finally, Fig. 8 shows an example of the behavior of the transverse rho propagator G p t- 
The mixing with the u meson is responsible for a moderate modification of the real and 
imaginary parts. 

6 Conclusion 

The dispersion relations of the neutral mesons a, u, p and 5 have been obtained in the random 
phase approximation in the framework of a relativistic hadronic field theory. It was shown 
that in nuclear matter with a proton-neutron asymmetry, new mixing channels open between 
mesons of different isospin number. A numerical evaluation indicates that the influence of 
this mixing on the effective meson masses and zero sound mode is moderate, of the order of 5 
% to 10 %. An application of these to the calculation of RPA corrections to neutrino-nucleon 
scattering also leads to effects of this order of magnitude. A similar study concerning the 
charged mesons is presently underway. 

The amount of corrections arising from meson mixing is somewhat larger in the p-uj 
(both transverse and longitudinal) and 5-u channels than in the p-a or a-5 channels for 
the thermodynamical and kinematic conditions explored in this work. The role of these 
new mixing channels on the behavior of the dilepton production rate in heavy ion collisions 
[§, [| remains to be investigated, and could bring further arguments to the discussion of the 
behavior of the p meson in the medium. 
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Appendix A - Real, Matter Part of Polarizations 



The general expressions of the matter part of the polarizations can be written in terms 
of the integrals for i £ {p, n} 



I 



m(i)(P,U,k) 



%k(i)(fo,k) 



du- 



(pu) n 



(cvjp 2 + Mf - pku - q 2 /2){toJp 2 + Mf - pku + q 2 /2 



p 2 dp 



X, 



(0 



p 2 + M 2 
p 2 dp 
p 2 + M 2 



(p 2 + M 2 ) n I m {p,ui, k) (rnlp) + (-l) n+m % 
(ni(p) +ni{pYj 



(53) 

(54) 
(55) 



Expressions for II CT , n^", 11^, n^, Wjf and 11^ have already been given in Appendix A of 
H in the case of symmetric nuclear matter. In asymmetric matter, one should remove the 
degeneracy factor d = 2 of the formulae presented in [|| and sum instead on contributions of 
the proton and neutron loops. For example, the S-p mixing polarization is now given by 



Tie U Sp 



TZeU 
95 



Sp(p) 



+ 1Ze n 



5p(n) 



f, 



P_ n "2\ n 2n 



^(2M i5p + ^V)^ 



% 



T- 







o(0 x iW 



(56) 



In asymmetric matter, we will moreover have mixing polarizations between mesons of different 
isospin number. They are given by the difference between the contributions of the proton 
and neutron loops. 



Tie UupL 



K-e U wp l ( p ) - Tie HujpL (n) 



TZeU 



uipL (i) 



9m9p 

' 2tt 2 



X, 



0(t) "^2(0 



A«\ - {^A ^M t q % 



\2m) 4tt 2 ' 



'0(0 



Tie UujpT = Tie U upT ( p ) - Tie Ti upT ( n ) 



Tie ~R U pT(i) 

TZeW ap 



erp{i) 

Tie U^X 



9u9 P 

~4vr 2 



Mf qX(i) + (^ + k 2 )lL) +I° 2(i) ~ 4ukXl 



Suj(p) 



'(0 

^ n L(n) 



!(0 



(0 



989<, 
' 2ir 2 



■if Mi q 2 



T 1 

H{i) 



w T o 



^ K( P ) 



^ Kin) 



-g a I 2Mi g p + ^-q 2 



2m 



4vr 2 



% 



o(0 



w T o 
t 2 i(0 



TZeU 



a5(p) 



TZeU 



TZeU 



<tS (i) 



9a98 

' 2n 2 



(0 



aS(n) 
A 



2XX> + 



M?^ (ll 
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Appendix B - Renormalization 

Let us first consider the scalar mesons: 

The contribution of the vacuum to the a meson polarization is 



Tie IC C 



flg 



f ) +r T - M; q 



with I 



(0 



X 



(0 



2tt 3 
d A p 



> ^i^+i^ + [\-MiAA n) 



, 4 6(p 2 -Mf)e(- Po ) 



(p.k) 2 - /fc 4 /4 



The vacuum part of the 5 polarization TZeHJ^ has the same structure. The divergent part 
of the integrals X\ and X\ is extracted by dimensional regularization 






(i) reg 



7rM 



r (i) reg 



2tt 



where 6 is a finite valued function 



9(q 2 ,M 2 )=y 



1 , /Afc . 

— + ln — Ai 
e \m 

1 , / M i A 
- + ln — A 2 
e V m 



cfa 



1 



+ 9(q 2 ,M 2 



with y = 1 



/o (x 2 +y)\/?TT ''' ' 4Mf 

In order to simplify somewhat the expressions, we will use the following short notations 



(57) 



6(q 2 ,M 2 ) , e n = 6(q 2 ,M 2 ) 
, 8s = 6(m 2 ,m 2 



e(mi,m 2 ) 

- 2 M 2 ) 



>oq 



dq 2 



\q 2 =in 2 ,M 2 =m 2 



(58) 



and work in units of the free nucleon mass (i.e., put m = 1). e is a vanishingly small quantity 
and Ai, A2 are finite arbitrary integration constants. 
We obtain infinities of the type 



IT 



2vr 2 



M 2 + M 2 q i 

—2 + — + A 

ei e 2 



with M 2 + M 2 = 2(m — g (J (y') 2 + 2g 2 5 2 and A CTCr is a finite contribution. Since the mixed 
polarization II^ C depends on the difference of the contributions of proton and neutron loops, 
one finds that the q 2 /e singularities from both contributions cancel each other, and one is left 
with 



nvac 
a5 



9u9S 
2vr 2 



M 2 
p 



£3 



M 2 



with M? — M 2 = 2(rn — g a a)gs5 and A a $ is a finite contribution. IIJ| C and II^ C have similar 
structures. 

These infinite contributions can be cancelled by adding to the Lagrangian (|3|) the coun- 
terterms 

Zeduvd^a + Zsd^S + \a 1( j 2 + l-A 2 a 3 + ]-A 3 a 4 + B x 5 2 + B 2 ab 2 + B 3 a 2 5 2 + C5 A (59) 
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which contribute to the dispersion relation as e.g. 

„2 



(-^ + M>! 



rrmat _|_ "(T 

aa + 2vr 2 



A, 



01 + 



■rrmat i 9cr98 



A n 



+ 



+ 



'zrf + A l + 2A 2 a + 3^ 2 + 2^ + A { *™^¥±Ml + £ ' 



2vr 2 



ei 



25 2 5 + 4S 3 o-f5 + 



ff CT fl<5 f 2(m - g a a)g s 5 
2tt 2 1 e 3 



<Ji 



The infinities are neutralized by choosing the couplings of the counterterm Lagrangian so 
that, for example 



o\ 



n 2 n 2 1 



45, 



2B 3 + 



vr 2 
9 2 A 


ei 
1 


TT 2 

gWs 


£3 

1 


IT 2 


£4 
1 



L<T<5 



45, - ^±— = bl° 

it 2 e 6 



zi>2 + m — ^ — On 

it 2 e 3 

-k 2 e 4 

9 R , __ 9a95 ! _ . <5fj 

zr>2 + <7i — 5 — o 2 

vr 2 e 6 



(60) 



with &3 CT , 6f , ••• etc finite. The unknown finite remnants are defined so that the vacuum 
polarizations 



-rrvac 

(JO 



nvac 



^Ka + z a q 2 + ai + 2a 2 a + 3a 3 a 2 + 2b% <7 5 2 

Z7T Z 

" "" A ffS + 2b? 5 + 46f 



2vr 2 



fulfill some physical conditions. One generally requires the polarization Il™ c to cancel in the 
vacuum (M p = M n = m, a = 5 = 0) on the mass shell q 2 = m 2 . In J28|] one further imposes 
that the first derivatives with respect to the a and 8 fields vanish at the same point 



d ff 



-n d2U ™ \ 

\q 2 =m%,M p =M n =m,cr=5=0 ~ u i n 2 \q 2 =m 2 r ,M p =M n =m,a=8=Q 







(61) 



In the scheme of Kurasawa and Suzuki (scheme A in this work), one requires that these 
derivatives vanish at q 2 = 



dm 



Oa 



\q2=0,M p =Mn=m,cr=d=0 







d 2 Ul 



"~n2 \q 2 =0,Mp=M n =m,<j=S=0 







(62) 



In the schemes of Shiomi and Hatsuda [33] or scheme 3 of |U| (scheme B of this work), the 
structure of the expressions in terms of M n , M p is preserved, leading to conditions such as 



01 + 2a 2 a + 3a 3 a 2 + 2b% a 5 2 = ai— ^ [2(m - g^af + 2g 2 s 5 2 

2m z L 



a\ 



2m 2 



M 2 + M 2 



(63) 



These conditions yield expressions of the constants 63 CT , 6f , ••• e ^ c m terms of the function 
9 defined in Eqs. (|5H). For example, in scheme B one obtains 
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The constants determined by these conditions must moreover fulfill the compatibility condi- 
tions imposed by (BH) 



bf-b S / = ^(b s /-bf) , 2b s 3 s -bf = ^(b?-b 5 2 5 ) (64) 

LCrS V&U L(5o 



The constants 63 , 63 , b% as determined by the sets of conditions (p2j) or (p3j) do indeed 
fulfill these compatibility conditions, whereas the set of conditions (|6l|) do not. 

We now give the full expressions of the renormalized polarizations in schemes A or B. 



Renormalization scheme A 



n 2 
j-rvac (A) _ "a 

2vr 2 



3M 2 lnM p + 3M 2 lnM n - ^ (lnM p + 9 P )-\ QnM n + n ) 

+2 M 2 9 p + 2 M^„ + [q 2 - ml) (o a - (4 - m 2 ,)^) - (4 - m 2 )0, 

/ M 2 + M 2 \ " 

+ 1 _ 13 " £ + 6(M n + M p ) (65) 



The expression of IlTt c can be obtained by replacing everywhere g a , m a 9 a , 9 aq by g$, m$ 



), "dq 



. -pj-vac (A) 

in n CTCT 



n 



vac (A) 



9ct98 

2tt 2 



3M 2 InMp -3M n 2 lnM n -| (lnM p + 9 P ) + y (lnM n + r> 



+2M p \ -2M„X + ( Mn n Mp ) (-12 + 13(M n + M p )) 



(66) 



Again, the expression of IT is obtained by replacing everywhere the indices a by 5 in 

TT VaC (A) 

For the vector mesons, we have 

J2 



■prvac (A) fiu 



9Z_ 

6tt 2 



2M 2 (9 p - 1) + 2M 2 (9 n -l)+q 2 (In M p + 9 p ) + q 2 (In M n + r , 



4g 2 



'"- 



(0 W - 1) - 2</ 2 



x \g 



AiU 



qHqf 



(67) 



■prvac (A) jxv 



■prvac (A) jiv 

2M 2 (9 P - 1) - 2M 2 (9 n - 1) + q 2 (In M p + p ) - q 2 (+ In M n + 9 n ) 



9u9 P 
6vr 2 



+ ( {£) 2^ q2 [Mp (ln Mp + 0p) " M " (ln M n + n ) + 2(M n - M p )] 



x 3- 



7 /i„i/ 



/W 



(68) 



-prvac (A) /iv 
PP 



6vr 2 



2M, 



p\ w p 



1) + 2M 2 (0 n - 1) + q 2 (In M p + 9 p ) + c? 2 (In M n + n ) 
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-^(9 p -l)-2q 2 



2 9 



+ 



(AY a 

\2mJ 6ir 2 

,2 



3M 2 In M p + 4M 2 fl p + |- (In M p + p ) + 3M 2 In M„ + 4M 2 0, 



+^- (lnM n + n ) - (8 + m 2 )0 p - (<? 2 - m 2 p )(6 p + (8 + m^g) 

M 2 + M 2 
+5 + 6(M n + M p ) - 17—2— 



+ 



/p \ 5p 9 



2m/ 7T 2 



M p (In Mp + p ) + M n (In M n + 9 n ) - 2d p + 4 ( m 

(1 i/> 



M n + M p 



nv _<rq_ 



(69) 



Renormalization scheme B 



-Qvac (B) 



2vr 2 



p 



(In M n + 9 r , 



:J M 2 lnM p + 3M 2 lnM n - % (lnM p + 6 P ) - 

2 z, 

+2 M 2 # p + 2 MX + 1 - p " (40 CT - m 2 (4 - m 2 )0„ 



+ (q 2 - ml) (6 a - (4 - m 2 )^) - (4 - m, 



1)0, 



(70) 



n 



vac (B) 



9a9S 

2vr 2 



3M p 2 kiM p - 3M 2 lnM n - \ (lnM p + 9 p ) + \ (lnM n + 0„ 



/M 2 - M 2 \ / 
+2M p X - 2 M 2 # n + p - (-40 CT + m 2 (4 - m 2 ) 



'<rq 



(71) 
(72) 



Once more, the expressions of n^ c and n^ c are obtained by replacing everywhere the 

indices a by 8 in 14^ and II^ C respectively. 

The expression of H^P" coincides in both schemes. For the polarizations involving a 
vertex with the rho meson, we have 



T-rvac(B) \iv 



9co9p 

6vr 2 



2M 2 (0 p - 1) - 2M 2 (0 n - 1) + q 2 (lnM p + 9 p ) - q 2 (+lnM n + 6 n ) 



+ ( ^£j T^q 2 [M p (lnMp + 6 p ) - M n (lnM n + 6 n ) + (M n - M p )9 w ] 



x \9' 



(73) 



and Ilp^f is obtained by replacing 9 W by P in the last term. 
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■prvac (B) fiu _ I 9p 

pp I 6vr 2 



2M 2 p (d p -l) + 2M*( 



1) + q 2 (In M p + e p ) + q 2 (In M n + n ) 



49 M 

-3 s , 



2 V V P 



+ 



m 



2 9 



2g 2 P 



2m / 67r 2 

2 



3M p 2 In M p + AM 2 9 P + ^- (In M p + 9 P ) + 3M 2 In M n + 4,M 2 9 n 



+ ^- (lnM n + n ) - (8 + m 2 )^ - {q 2 - m 2 )(6 p + (8 + m 2 )6 pq ) 



+ 



' M 2 + M 2 ^ 



+ m 2 p (8 + m 2 )0 pg 



+ 



^-) ^q 2 [M p (In Mp + 9 p ) + M„ (In M n + O n ) - 9 p (M n + M p )) 



x \9 



flU 



lf i n u 



(74) 



As mentioned before, the renormalization constants fulfill the compatibility conditions 
(|64[). Note however that, in scheme B, H a $ 7^ n^ since they are renormalized on the mass 
shell of the a and 5 respectively, i.e. 



lUte 2 



m~ 



;n 5ff (g 2 



m 2 ) 







This breaks the symmetry which was explicit on the initial divergent expression. The symme- 
try could be reestablished e.g. by imposing that the polarizations vanish at some intermediate 
common scale fi a $ 

Ras{q 2 = vl&) = = H-sa(q 2 = /4<s) 

In renormalization scheme A on the other hand the symmetry was already preserved by the 
choice of the common renormalization point q 2 = 0. The same remark applies to 11^ 7^ IP^ 
in scheme B since they are renormalized on the mass shell of the 00 and p respectively. They 
coincide in scheme A. 



Appendix C - Imaginary Parts of Polarizations 



The imaginary part arise the factor ±ie in the denominator of ([24]) after applying the 
formula l/(x + ie) = V jx — itt5(x). Depending on the prescription applied for going around 
the pole in the propagator, we obtain the imaginary part of one of the definitions n, U R , n < , 
II 11 ... all these being related to each other through factors such as tanh(/3u;/2) pG|1 . Here 
we chose to give the retarded polarizations. They are obtained by inserting ±ie sign(po) in 
the denominator of Eq. ( p4[) . 

All imaginary parts can be expressed in terms of three integrals E[ , E% , E$ with 
i £ {n,p}. At finite temperature, for the calculation of the retarded polarizations, we have 
— for spacelike momentum: 



E (i) 



dy 



V+ 2 



n-l 



+(-i) n y 



o(y - vl ){ n i(.y) - My + w )} 



, > re— 1 ,.. 

% - Vu) {ni{y) - ni(y - u>)} 



(75) 
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— for timelike momentum with q 2 < 4M 2 , the imaginary parts vanish: E n = 

— for timelike momentum with q 2 > AM 2 : 



E (i) 
J ~'n 

with 



\n-l 



UJ 



n—1 



.«> 



u dy [ y ~ I j i 0(y " ^ " e{y ~ $)] { n(i) ^ - y) + ™ {i) (y) -(* 



niiy) 



Vl 



kVA® - to 



-i 



«i(y) 



Mv+Vi) _i_ i 



-i 






fcVAW +uj 



A (i) = 1-4 



Af? 



Xm nW 






Imff 



W 



Xm n 



Ira nW 



2vr/c 

g 2 

47r/c 



^ - ^ E 



(0 



n« 



<? 2 n« 



Mf + V | B^ + ^£ 3 



fc 2 " 



Xm n 



(0 



"2^3 1 M * ^ 

._?L (m 2 - ^-) E 

2nk \ l 4 



Ira n 



(0 

pL 



1 



Xm n 






2nk 
1 

47T/C 



,i(0 



W 

1 

2 



T E i ' - F E 



(0 



+ ,^ g =»M i ^ + (A N2 



'2m 



2m 



M 2 g 2 B « + |_^) 



9 2 p 



^+Tl^+^ 



(i) _l <T p(i) 
/fc 2 



+ 2 ff A g 2 M . £,(») 



2m 



+ 



2m 



r2„2 , T 



n(i) 



9 4 ™(0 



^ + ^ #r - ^ 



Xm n 



(0 

5p 



Xmng 



gg (2Mj gp + {f P /2m)q 2 ) q 2 (i) 

4vrF ^ 2 

gCT (2M f g p + (/ P /2m)g 2 ) g 2 (i) 



Xm, n 



Xm n 



(0 

w 



ink 3 
9s9ui . , 2 P (») 

-M^f M 2 _^Uw 

2vryfc I * 4 / 1 



Eh 



Zm n 



(0 

pLU L 



1 



Xm n 



(0 

PT 



Airk 

1 

Airk 



9 P 9u 
9 P 9u 



-E\ 



(0 



4^ 



+9uj I^ q 2 M i E? 



W4)^4^ 



+ <^g 2 M,£« 
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Fig. 1 Effective masses and pro- 
ton fraction as calculated from pa- 
rameter set ( |49| ) with gs = 5 in 
neutrino-free matter in (3 equilib- 
rium at T=0. 
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Fig. 2 Effective masses of mesons as calculated with renormalization procedure A as a function of 

asymmetry. The density was fixed at n# = 4n sat and the temperature at T = 0, and the coupling 

constants are those of (fj 
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Fig. 4 p-w mixing angle in the transverse mode as a function of density asymmetry, as calculated 
with renormalization procedures A (left) or B (right). The parameters are those of [Q2J and the 
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Fig. 5 Effect of asymmetry on the 
zero sound mode due to a-Lu mixing. 
The figure was plotted for ub = 4n sat 
and T = and using renormalization 
scheme (A) . The coupling constants are 
those of |22|. 
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Fig. 8 Real and imaginary parts of the transverse p meson propagator. 
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